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Almost global existence for the Prandtl boundary layer equations 

Mihaela Ignatova and Vlad Vicol 

Abstract. We consider the Prandtl boundary layer equations on the half plane, with initial datum that lies 
in a weighted space with respect to the normal variable, and is real-analytic with respect to the tangential 
variable. The boundary trace of the horizontal Euler flow is taken to be a constant. We prove that if the Prandtl 
datum lies within e of a stable profile, then the unique solution of the Cauchy problem can be extended at least 
up to time Te > exp(e“^/ log(£“^)) . 


1. Introduction 


We consider the two dimensional Prandtl boundary layer equations for the velocity field , v^) 

dtu^ — dyU^ + u^dxU^ + v^dyU^ = —dxP^ ( 1 . 1 ) 

dxU^ + dyV^ = 0 ( 1 . 2 ) 

posed in the upper half plane BI = {(x, y) G : y >0}. Here denotes the trace at clIHI of the underlying 
Euler pressure. The boundary conditions 

U^\y=Q = V^\y=Q = 0 (1.3) 

U^\y=oo = (1-4) 

are obtained by matching the Navier-Stokes no-slip boundary condition = 0 on clH, with the Euler slip 
boundary condition at y = oo. The trace at OH of the Euler tangential velocity u®, obeys Bernoulli’s law 

dtu^ + u^dxu^ + dxP^ = 0 . 


The Prandtl system (1.1)-(1.4) is supplemented with a compatible initial condition 

u^\t=o = Uq. (1.5) 

Our main result states that if the Euler data {u^ ,p^) is constant, and if the initial datum of the Prandtl 
equations lies within e of the error function erf(y/2) (in a suitable topology), then the Prandtl equations 
have a unique (classical in x weak in y) solution on [0, Tg], where Tg > exp(e“^/ log(e“^)). 

Theorem 1.1 (Almost global existence). Let the Euler data be given by = k and dxP^ = 0. 
Define 

uo{x, y) = Uq {x, y)- k erf (JQ 

where erf is the Gauss error function. There exists a sufficiently large universal constant C* > 0 and a 
sufficiently small universal constant e* > 0 such that the following holds. For any given e G (0, e*], assume 
that there exists an analyticity radius tq > 0 such that 


logi 


<3/2 


< 


1 


1 
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and such that the function 

9oix,y) = dyUo{x,y) + ^uo{x,y) 

obeys 

II9oIU2,^,i/ 2 := ||exp(^)ag^go(x,i/)||^2(H)(2ro)"^ ^^^ - < e. 

m>0 

Then there exists a unique solution of the Prandtl boundary layer equations on [0, Tg], where 

The solution is real analytic in x, with analyticity radius larger than tq /2, and lies in a weighted 
space with respect to y. We emphasize that e and tq are independent of k. 

The precise function spaces, in which the solution lies, are given in Theorem 2.2 below. The condi¬ 
tion relating s and tq stated above roughly speaking says that we think of 0 < e ^ 1, and of tq = 0(1). 
The stated condition is the sharp version of this heuristic. 

Remark 1.2 (Initial vorticity may change sign). We note that the initial datum Uq is not necessarily 
monotonic in y, i.e. we do not necessarily have loq := dyU^ > 0 or < 0 on H. Thus, the initial data 
in Theorem 1.1 need not fit in the Oleinik [01e66] sign-definite vorticity setting. To see this, one may for 
example consider k = e > 0 sufficiently small and tq = 1/4. We then let 

«o {x, y) = e (exp(-x2)r/(?/) -|- erf(i//2)) , 

with p(y) such that p(0) = 0 and ex.p{y‘^ / A)q(y) G Then 

dyUo(0,y) =£(p'(y) -hexp(-yV4)/V^) 

can be designed so that 

dyUQ(0,0) > 0 and dyUQ(0,l) < 0- 

This indeed shows that the initial profiles considered in Theorem 1.1 need nof be monofonic in y. 

1.1. The local well-posedness of the Prandtl equations. Before discussing the proof of our main 
result (cf. Subsection 1.3 below), we present the history of the problem. The Prandtl equations arise from 
matched asymptotic expansions [Pra04] meant to describe the boundary behavior of solutions to the Navier- 
Stokes equations with Dirichlet boundary conditions 

+Vp^’^'> = 0, inn, (1.6) 

= 0, on dn 

in the vanishing viscosity limit —)• 0. Here C is a smooth domain. Formally, as z/ —)• 0 the Navier- 
Stokes equations reduce to the Euler equations, for which the slip boundary condition ■ n = 0 holds on 
dn. Due to this mismatch of boundary conditions, uniform in u bounds for in e.g. the L^(n) norm, 

on an 0(1) time interval, remain an outstanding mathematical challenge. 

One of the fundamental questions which arise is to either prove that the Prandtl asymptotic expansion 

(u^''\v^''^)(t, X, y) = (u^, v^)(t, X, y) + (u^ , y/uv^){t, x, y/s/u) + o(v^), (1.7) 

can be justified rigorously [SC98, Mael4, GN14], or to show that it fails [GreOOa, GVDIO, GNU, GVN12, 
GGN14b, GGN14c, GGN14a]. Naturally, the answer is expected to depend on the topology in which (1.7) 
is considered, and this is intimately related to the question of well-posedness of the Prandtl system. By 
now, the local in time well-posedness of (1.1)-(1.5) has been considered by many authors, see e.g. [01e66, 
EE97, SC98, CSOO, GreOOb, CLSOl, HH03, LCS03, GSS09, MW14, KV13, GVM13, Mael4, AWXY14, 
KMVW14, LWY14, WXY14, LWX15] and references therein. However, the question of whether the 
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inviscid limit holds whenever the Prandtl equations are locally well-posed, and are thus stable 

in some sense, remains open. See [CKV14] for partial progress in this direction. 

In [01e66], Oleinik proved the existence of solutions for the unsteady Prandtl system provided the 
prescribed horizontal velocities are positive and monotonic, i.e. > 0 and lo^ = dyU^ > 0. From the 
physical point of view, the monotonicity assumption has stabilizing effect since it prevents boundary layer 
separation. The main ingredient of the proof is the Crocco transform which uses as an independent 
variable instead of y and as an unknown instead of . More recently in [MW14], Masmoudi and Wong 
use solely energy methods and a new change of variables to prove the local in time existence and uniqueness 
in weighted Sobolev spaces, under the Oleinik’s monotonicity assumption. The main idea of [MW14] is 
to use a Sobolev energy in terms of the good unknown — u^dy log(w^), which may be done if 

u;P>0. The equation obeyed by the top derivative in x of is better behaved than that of the top derivative 
of either or Although cf. Remark 1.2 this change of variables is unavailable to us, the idea of a good 
unknown inspired by [MW14], plays a fundamental role in our proof. Recently, in [KMVW14], the local 
existence and uniqueness for the Prandtl system was proven for initial data with multiple monotonicity 
regions, as long as on the complement of these regions the initial datum is tangentially real-analytic. 

For real-analytic initial datum, Sammartino and Caflisch [SC98, CSOO] established the local well- 
posedness by using the abstract Cauchy-Kowalewski theorem. For initial datum analytic only with respect 
to the tangential variable the local well-posedness was obtained in [LCS03]. In [KV13], the authors gave an 
energy-based proof of this fact, and considered initial data with polynomial rather than exponential match¬ 
ing at the top of the boundary layer u^{t, x, y) — u^{t, x) —)• 0 as oo. The tangentially-analytic norms 
introduced in [KV13] encode at level a full one-derivative gain from the decaying analyticity radius. 
These norms play an essential role in our proof. 

More recently, for data in the Gevrey class-7/4 in the tangential variable, which has a single curve of 
non-degenerate critical points (i.e. ujq = 0 iff y = ao(x) > 0 with dyUjQ (x, ao(x)) > 0 for all x), the local 
well-posedness of the Prandtl equations was proven by Gerard-Varet and Masmoudi in [GVM13]. Note that 
this Gevrey-exponent is not in contradiction with the ill-posedness in Sobloev spaces established by Gerard- 
Varet and Dormy in [GVDIO] for the linearized Prandtl equation around a non-monotonic shear flow. Here 
the authors show that some perturbations with high tangential frequency, k ^ 1, grow in time as At 
the nonlinear level this strong ill-posedness was obtained by Gerard-Varet and Nguyen [GVN12]. 

1.2. The long time behavior of the Prandtl equations. As pointed out by Grenier, Guo, and Nguyen 
[GGN14a, GGN14b, GGN14c] (see also [DR04]), in order to make progress towards proving or disproving 
the inviscid limit of the Navier-Stokes equations, a finer understanding of the Prandtl equations is required, 
and in particular one must understand its behavior on a longer time interval than the one which causes the 
instability used to prove ill-posedness. However, to the best of our knowledge the long-time existence of the 
Prandtl equations has only been considered in [01e66], [XZ04], and [ZZ14]. 

Oleinik shows in [01e66] that global regular solutions exist, when the horizontal variable x belongs to a 
finite interval [0, L], with L sufficiently small. Xin and Zhang prove in [XZ04] that if the pressure gradient 
has a favorable sign, that is dxP^{t, x) < 0 for all f > 0 and x E M, and the initial condition Uq of Prandtl 
is monotone in y, the solutions are global. However these are weak solutions in Crocco variables, which 
are not known to be unique or to be regular. The global existence of smooth solutions in the monotonic 
case remains to date open. In the case of large datum, the assumption of a monotone initial velocity is 
essential. Indeed, E and Engquist [EE97] take = dxP^ = 0 and construct an initial datum Uq which 
is real-analytic in the tangential variable, but for which ojq is not sign definite, and prove that the resulting 
solution of Prandtl (known to exist for short time in view of [LCS03, KV13]) blows up in finite time. We 
emphasize that for this blowup to occur, the initial datum must be at least 0(1): indeed, for initial datum 
that is sufficiently small, the conditions of Lemma 2.1 in [EE97] fail, and thus the proof does not apply. 

In fact, for initial datum that is tangentially real-analytic and small, in a recent paper Zhang and 
Zhang [ZZ14] prove that the system (1.1)-(1.5) has a unique solution on a time interval that is much longer 
than the one guaranteed by the local existence theory. More precisely, for = e and dxP^ = 0 it is proven 
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in [ZZ14] that if Uq = 0{e) in a norm that encodes Gaussian decay as ?/ — )• oo and tangential analyticity 
in X, and if e <C 1, the time of existence of the resulting solution is at least The elegant proof 

relies on anisotropic Littlewood-Paley energy estimates in tangentially analytic norms, inspired by the ones 
previously used by Chemin, Gallagher, and Paicu [CGPll] to treat the Navier-Stokes equations with datum 
highly oscillating in one direction (see also [PZ14] and references therein). 

1.3. Almost global existence for the Prandtl equations. In [ZZ14, Remark 1.1], the authors raise the 
following question: “whether the lifespan obtained in Theorem 1.1 is sharp is a very interesting question”. 
That is, do the solutions of the Prandtl equations with size e initial datum live for a time interval longer 
than 0(6“^/^)? In this paper we give a positive answer to this question, and prove (cf. Theorem 1.1 or 
Theorem 2.2) that in 2D we have almost global existence (in the sense of [Kla83]). That is, the solution 
lives up to time C)(exp(e“^/log e“^)). Our initial datum Uq consists of a stable 0{k) boundary layer lift 
profile, and an 0{e) possibly unstable, but tangentially real-analytic profile. In particular, the total initial 
vorticity is not necessarily positive (cf. Remark (1.2)). Whether solutions arising from sufficiently small 
initial datum are in fact global in time remains open, and this may depend on whether k < e or e <C k. 

The proof of Theorem 1.1 proceeds in several steps. In order to homogenize the boundary condition at 
y = oo, we write as a perturbation tt of a stable shear profile y), with dy(p{t, y) > 0. In order to 
capture the maximal time decay from the heat equation, and to explore certain cancellations in the nonlinear 
terms of Prandtl, we choose the boundary lift ip{t, y) to be the Gauss error function eY{{y/y/ 4(f + 1)). The 
equation obeyed by u (cf. (2.5)) contains the usual terms in the Prandtl equations, but also terms that are 
linear and quadratic in ip. The lift ip is chosen so that the quadratic terms in ip vanish, and we are left to 
understand the linear ones in ip. We note that until here a similar path is followed in [ZZ14], and that energy 
estimates for the ensuing linear problem leads to the maximal time of existence 0(e“^/^). 

The main enemy to obtaining a longer time of existence is the term Kvdyip in the velocity equation 
(2.5) (respectively Kvdyip in the vorticity equation (2.7)). As v = —dy^dxU, this term loses one tangential 
derivative, and is merely linear in u, so that it is not small with respect to e. 

The main idea of our paper is to introduce a new linearly-good unknown g = lo — udy{logdyip), 
cf. (2.12) below. This change of variable is directly motivated by the one in [MW12] for the case ui^ > 0. 
The upshot is that g obeys an equation in which the bad terms nvdyip and nvdyip cancel out, cf. (2.20) 
below. The solution of this new equation may be shown to be globally well-posed. Note here that we may 
recover u = \]{g) and v = \{g) via linear operators U and V that are nonlocal in y, cf. (2.23)-(2.24). Thus 
g is the only prognostic variable in the problem, and the system (2.20)-(2.22) is equivalent to (1.1)-(1.4). 

In order to take advantage of the time decay in the heat equation, it is natural to replace the y variable 
with the heat self-similar variable z = z{t, y) = y/y/t + 1, and to use norms with gaussian weights in 
the normal direction. The gaussian weights are useful when bounding u and v in terms of g, cf. Lemma 3.1. 
Moreover, the gaussian weights allow us to deal with another technical obstacle; namely, that in unbounded 
domains the Poincare inequality does not hold. However, with the gaussian weights defined in (2.25), 
we may use a special case of the Treves inequality (cf. Lemma (3.3)) as a replacement of the Poincare 
inequality. The need for a Poincare-type inequality in the y variable comes from the desire to work with 
norms, and still capture the full one-derivative gain from the decay of the analyticity radius. As was shown 
in [OTOl, KVll, KV13] this may be achieved by designing norms based on rather than sums over the 
derivatives. The one derivative gain inherent in these -based norms allows for direct energy estimates. 

The main ingredient of the proof is the a priori estimate (3.23) below. The idea is to solve the PDE 
(2.20)-(2.22) for g simultaneously with a nonlinear ODE (3.24) for the tangential analyticity radius r. The 
fact that the analyticity radius r does not decrease to less than r(0)/2 on the time interval considered follows 
from the time integrability of the dissipative terms present on the left side of (3.23) (see also [PVll]). 

The proof of Theorem 1.1 is concluded once we establish the uniqueness of solutions in this class, and 
show that there exists at least one solution to the coupled system for g and r. While uniqueness follows from 
the available a priori estimates, the existence of solutions introduces a number of additional difficulties. One 
of these is proving existence of solutions to (3.24). This is a first order ODE in r, for which the nonlinear 
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forcing term is well defined (i.e., the infinite sum converges), only if the solution g already is known to have 
analyticity radius r. To overcome this difficulty, we consider a dissipative approximation of (2.20)-(2.22) 
and for u > 0 add a —vd^g term on the left side of (2.20). We prove that this regularized equation has 
solutions g^'^'> in a fixed order Sobolev space, and a posteriori show that these solutions are tangentially 
real-analytic, with radii that obey an ODE similar to (3.24). Here we essentially use that the initial 
datum po is assumed to have tangential analyticity radius 2ro, while the solution is only shown to have a 
radius that lies between tq/2 and tq. We prove that these radii are uniformly equicontinuous in v (in 
fact uniformly Holder 1/2) so that they converge along a subsequence on the compact time interval [0,Te]. 
To conclude the proof, we show that along this subsequence the are a Cauchy sequence when measured 
in the tangentially analytic norms, and that the limiting solution g and limiting radius r obey (3.24). 

Organization of the paper. The detailed reformulation of the Prandtl system is given in Section 2. 
Here we also define the spaces in which the solutions lives, and reformulate Theorem 1.1 in these terms. 
The a priori estimates are given in Section 3, the uniqueness of solutions is proven in Section 4, and the 
details concerning the existence of solutions are given in Section 5. 


2. The linearly-good unknown, function spaces, and the main result 


Denote by 


z = z{t,y) = 


where (t) = t + 1 


the heat self-similar variable. We consider the lift 

Kif = Kip{t,y) 

of the boundary conditions (1.3)-(1.4), where 

ip{t,y) = 4>(z(t,y)) 

and the function 4>(2;) obeys 

4>(0) = 0, lim < 1 >( 2 ;) = 1, ^'{z) > 0. 


(2.1) 


( 2 . 2 ) 

(2.3) 


We make a precise choice of in (2.18) below. We already note that by design dy(p{t, y) > 0 for y > 0, 
i.e., the vorticity of the shear flow is positive, and thus stable (in the sense of [01e66]). 

We write the solution of (1.1)-(1.4) as a perturbation u{t, x, y) of the lift Kip{t, y) via 

X, y) = Kip{t, y) + n(t, X, y) 

so that the perturbation u obeys the homogenous boundary conditions 


— '^\y=oo — fl? 

and satisfies the equation 

dtu — dyU + KifdxU + Kvdyip + udxU + vdyU = — dy^p), 

where v is computed from u as 

ry 

v{t,x,y) = - dxu{t,x,y)dy. 

Jo 

Using (2.5), we obtain the equation for the perturbed vorticity uj = dyU, 

dtU! — dyOj + KifdxUj + Kvdyip -|- udxOJ + vdyU = —K{dtdy(p — dy(p), 

with the natural boundary condition 


(2.4) 

(2.5) 

( 2 . 6 ) 


(2.7) 


dyUj\y=0 = -Kdy(p\y=0. 


(2.8) 
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As we work in weighted spaces, at y = oo we impose the condition 


(2.9) 


2.1. The linearly-good unknown. As can already be seen in [ZZ14], the main obstruction for obtain¬ 
ing the global in time existence of solutions comes from the linear problem for the velocity and vorticity 

dtu — dyU + KifdxU + Kvdyip = —K{dt^p — dy^) + nonlinearity, (2.10) 

dtU! — dyOJ + KifdxOJ + Kvdyip = —K{dtdy(p — dy(p) + nonlinearity. (2.11) 

Inspired by [MW12] (see also [GVM13, KMVW14]), we tackle this issue by considering the linearly-good 
unknown 

9{t, X, y) = u{t, X, y) - u{t, x, y)a{t, y) ( 2 . 12 ) 


where 


a{t,y) = 


dyV>{t,y) 


dyip{t, y) ■ 

Note that one may solve the first order (in y) linear equation (2.12) to compute u from g explicitly as 

fy _ I _ n _ I _ 

uit, X, y) = dyip{t, y) / g{t, x, y)-^ — TT^^y = 2/)) / 9(t, x, y) dy, 

Jo Oy(p{t,y) Jq ^’{z{t,y)) 


(2.13) 


(2.14) 


where we have used the boundary condition of u at y = 0. Also, if g decays sufficiently fast at infinity, 
this ensures the correct boundary conditions for u. The formula (2.14) is useful when performing weighted 
estimates for u in terms of weighted norms of g. 

The evolution equation obeyed by prognostic variable g is 


dtg — dyg + (n + rnp)dxg + vdyg — 2dyag + uL + vudyO = kF, 


(2.15) 


where the diagnostic variables u and v may be computed from g via (2.14) and (2.6). The functions F and 
L are given by 

F{t, y) = aidtip - d^cp) - (dtdyip - d^cp) 

L{t, y) = dta — dyO — 2adya, 


and a is as defined in (2.13). Moreover, in view of (2.4), (2.8), and (2.9), the linearly-good unknown obeys 
the boundary conditions 

{dyg -h ag)\y=o = dyUj\y=o = -Kdyip\y=o, (2.16) 


and 


9\y=oo — 0 , 

where the latter one comes from the convenience of vorticity that vanishes as y —oo. 


(2.17) 


2.2. A Gaussian lift of the boundary conditions. At this stage we make a choice for the boundary 
condition lift 4). Our choice is determined by trying to eliminate the forcing term F on the right side of 
(2.15), and the linear term uL on the left side of (2.15). For this purpose, let 4> be defined via 4>(0) = 0 and 


4>'(z) = —= exp 
/vr 



(2.18) 


where the normalization ensures that —>• 1 as 2 ^ oo. In the original variables, this means that 
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where erf is the Gauss error function. With this choice of <h, we immediately obtain 


a{t,y) = -- 


2(f) 


F{z) = L{z) = 0 

and the boundary values 

$"(0) = 0 and <h'(0) = > 0. 

The evolution equation (2.15) for the good unknown 


g = ui + 


2(f) 


u 


thus becomes 


1 


1 


dtg - dyQ + (n + Hip)d:cg + vdyg + —g - = 0 

dyg\y=Q = g\y= 0 O — 0 

g\t=o = go- 

As noted before (cf. (2.14) and (2.6)), u and v may be computed from g explicitly 


u{t,x,y) = \J{g){t,x,y) := exp - 


y 

4(f) 


g{t,x,y)exp 


f 

4(f) 


dy 


v{t,x,y) = V{g){t,x,y) := - \J{d^g){t,x,y)dy, 


(2.19) 


( 2 . 20 ) 

( 2 . 21 ) 

( 2 . 22 ) 


(2.23) 


(2.24) 


and thus solving (2.20)-(2.22) is equivalent to solving the Prandtl boundary layer equations (1.1)-(1.5). 

2.3. Tangentially analytic functions with Gaussian normal weights. Lastly, in view of (2.14) and 
the choice (2.18) of <h, it is natural to use the Gaussian weight defined by 


0a{t,y) = exp 


az{t,yf 


= exp 


ay 

4(f) 


(2.25) 


for some 


a G [1/4,1/2] 

to be chosen later (e-close to 1/2). 

In order to define fhe functional spaces in which fhe solution lies, mofivafed by [KV13], if is convenienf 
fo define 

+ 1 


Mm = 

and infroduce fhe Sobolev weighfed semi-norms 


m\ 


ig-iX) 


(2.26) 

Dm — Dm (y) t) 

1 = \\eadyd^g\\L2T'^Mm = Xm{dyg,T), 

(2.27) 

Zm — ^m{gj x') 

= \\z9ad'^g\\L2T'^Mm = Xm{zg,x), 

(2.28) 

Dm — Dm{g-i X ) 

= (f)^/^Xm(y,r) -h (f)^/^Zm(y,r) -h {tf^'^Dmig,x) 

(2.29) 

Tm — ^m{,g^x) 

= \\ead^g\\L2T'^~^mMm. 

(2.30) 


As in [KV13], we consider fhe following space of funcfion thaf are real-analytic in x and lie in a weighfed 
space wifh respecf fo y 

Xr,a = {g{t,x,y) G L‘^{B.;eadydx): ||y||x,,„ < oo} 
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where for r > 0 and a as above we define 

Ibik... = (2.31) 

m>0 

We also define fhe semi-norm 

\\9U^. = Y.^m{9,T) (2.32) 

m>l 

which encodes fhe one-derivafive gain in fhe analyfic esfimafes, when fhe summation in m is considered in 
rafher fhan in as is classical when using Fourier analysis. Nofe fhaf for /3 > 1, we have 

\\g\\Yr,c < r~^\\9\\xp^,^ sup (m/3"™) < CpT-^\\g\\xf,^^^. (2.33) 

m>l 

In particular, g € X 2 T,a implies fhaf ||( 7 ||y.,. „ < r"^ 11511 x 2 ^ a - The gain of a y derivative shall be encoded in 
fhe dissipafive semi-norm 

II^IId,,,. = ^^i9,r) = \\dyg\\x^,^, 

m>0 

while fhe damping in fhe heal self-similar variable z is measured via 

llffbr.c = ^rni9,T) = \\zg\\x^,^- 

m>0 

For compaclness of nofafion, for a function g such fhaf g, zg, dyg G Xr^a we use fhe lime-weighled norm 
I| 9 ||B.,. = E ( 2 . 34 ) 

m>0 


where as before r > 0 and a G [1/4,1/2]. Laslly, in order lo obfain time regularity for fhe radius of 
analylicily T(f), if will be convenienl fo use a hybrid of fhe 1“^ and fangentially analyfic norms, given by 


Iloilo.,, c 

:= X] 

m>0 

Dm -- 

= Dm{g,T) 

Dmid^rf 
Xmig,r) ’ 


(2.35) 


:= 

m>0 

Zm = 

- Zn 

ii9,T) 

_ Zm{9,Tf 
Xm{9,T) ’ 


(2.36) 


■= Bm 

Bm = 

= Bm{g,T) 

= {ty^'^Xmig: 

,t) + {t)^/'^Zmig,T) + {tf^‘^Dmig,T). 

(2.37) 


m>0 







We note that the bound 











WhWBr. 


1/2 II ||l/2 

(2.38) 


is an immediate consequence of fhe Cauchy-Schwarfz inequalify. 


2.4. The main result. Having inlroduced fhe funclional selling of Ibis paper we reslale Theorem 1 . 1 in 
Ihese terms. Firsf, we give a definition of solutions fo fhe reformulated Prandll equations (2.20)-(2.22). 


Definition 2.1 (Classical in x weak in y solutions). For /3 > 0 define to be the closure under 

the norm 

dydx 

m—n o—n 4HI \ / 


of the set of functions 


m=0 j=0 ' 

V = {h{x,y) G (^“(M X [0,(X))): dyh\y=o = 0}. 
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Let a G [1/4,1/2], and Oa(t, y) be defined by (2.25). For T > 0 we say that a function 

is a classical in x weak in y solution of the initial value problem for the Prandtl equations (2.20)-(2.22) on 
[0, T), if (2.20) holds when tested against elements of T) x M x [0, oo)). 

Theorem 2.2 (Main result). Assume the trace of the Euler flow is given by = k and dxP^ = 0. 
For t > 0, define 


u{t,x,y) = u^{t,x,y) — Kerf 


VW) 


y 


9{t, X, y) = dyu{t, X, y) + x, y). 


and let 

n(F nr* n i\ - / , , 

' 2(f) 

There exists a sufficiently large universal constant C* > 0 and a sufficiently small universal constant e* > 0 
such that the following holds. Assume that there exists an analyticity radius tq > 0 and an e & (0, e*] such 
that 


<T?P< ' 


(2.39) 


log 1 ° ’ 

and such that the initial condition go = g{0, ■, ■) is small, in the sense that 

119011x2,^,1/2 < e. 

Then there exists a unique classical in x weak in y solution g of the Prandtl boundary layer equations 
(2.20)-(2.22) on [0, Tf\ which is tangentially real-analytic, and the maximal time of existence obeys 


(2.40) 


Tg > exp 


--1 


log(e 


- 1 ' 


Moreover, letting 6 = elogl and a = the tangential analyticity radius r(f) of the solution g{t) 

satisfies 

2/3 


and the solution g{t) obeys the bounds 


X, 


i: 

L 


(*).< 




(s)V2 


> / 3/2 ' 

“T 21oglJ 

“ 2 

(2.41) 



(2.42) 

Bxis),a 4" 

log- 

(2.43) 

-||9(s)IIw(.).,II9(s)IIb,i 

, , ds < C^e 

(2.44) 


for all t G [0, T^]. 


It follows from the estimates in the next section (cf. Lemmas 3.1 and 3.5) that bounds on g, zg, and dyg 
in Xx a imply similar bounds on u and v in X, and thus (2.42)-(2.44) directly translate into bounds for 
uP and v^. Moreover, when g{t) G TL 2 ,i,al{t)’ thcri u{t) lies in T-L 2 , 2 ,al{t) the Prandtl equations (1.1) 
hold pointwise in x and in an sense in y. We omit these details. The proof of Theorem 2.2 consists of 
a priori estimates (cf. Section 3), the proof of uniqueness of solutions in this class (cf. Section 4), and the 
construction of solutions (cf. Section 5). 
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3. A priori estimates 

In this section we give the a priori estimates needed to prove Theorem 2.2. We start with a number of 
preliminary lemmas, which lead up to Subsection 3.5, where we conclude the a priori bounds. 

3.1. Bounding the diagnostic variables in terms of the prognostic one. We may use (2.14) to write 


0a{y)u{y) = f g{y)6a{y) 
Jo 


exp 


(1 a) ^ 2 „,2 


m 

On the one hand, it is immediate from the above that 

\\ 0 aU\\L'^ < PagWhl- 

On the other hand, for p G [1, 2] we may estimate 


{y -y)]dy. 


(3.1) 


(3.2) 


\0a{y)u{y)\ < \\ 0 ag\\LP/{v-i) 


y /p(l-a) 2 2 


exp 


4(f) 


{y -y )] dy 


i/p 




i/( 2 p) f DawsonF[ 2 :(f, y)Kp^a] \ 


{ 


K, 


p^a 


where Kp^a = v^p(l — a )/2 and 

ry ry 

DawsonF[i/] = exp(-y^) / exp(p^)dy = / exp(y^ - y‘^)dy. 

Jo Jo 


It is not hard to check that 


DawsonF[y] < 


1 + 2 / 

for all y 0. Because there exists a universal constant (JJ ^ 0 such that ly ^C + a + for p G [^1, 2j and 
a G [1/4,1/2], it follows that 

1 

\TP"~ " - 


\0a{t,y)u{t,x,y)\ < C{t)^/^‘^y">\\eag\\jp/(p-i)- 


(3.3) 


{l + z{t,y)Y/P 

for p G [1, 2]. Using (3.3) and recalling the definition of v in (2.6) we may prove the following estimates. 

Lemma 3.1 (Bounds for the diagnostic variables). Let 9a be given by (2.25) with a G [1/4,1/2], 
Define u = \J{g) and v = Y{g) by (2.23) respectively (2.24). For m>t)we have 


IlSoarsIlij < c(t)^i‘\\e^dT3 




)'J 


'ii,. 




\"F P^S^S,jg\td +C(tf^\\e^d"'g\\'ff \\ze„d’^gtff 

x,y t£,y tr,y tE,y 

\t P„3T*'g\\'d \\0^(fi‘d,gtd ||9„8r‘a»9lll/ 

x,y tr,y 

+ c{t)r^p„(t;g\\^iy kwrsllf.* Ik9a4"*+'9llff 

x^y tc,y x,y tc^y 

•^’fJ ■^x,y 

for some universal constant C > 0, which is independent of a G [1/4,1/2]. 




(3.4) 

(3.5) 

(3.6) 

(3.7) 

(3.8) 

(3.9) 
(3.10) 
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Proof of Lemma 3.1. From identity (3.1) we have 

\d^u{y)\ < \eo,{y)d^g{y)\exp 


(1 - a) 


4(t) 


(y -y )]dy 


<C{t)'/*\\e^(Cgh,. 

The bound (3.4) follows by taking the norm in x of the above, while the bound (3.5) follows upon 
additionally applying the ID Agmon inequality in the x variable, 

To bound Oad^u, we note that for i? > 0 we have 


IKd^ghi = 


rR 


IR 


\yeaiy)d^g{y)\\y\ ^dy 


\da{y)d^g{y)\dy + 

Jo 

< R^I^\\e^d^g\\Li+ R-^/^\\ye^d^g\\Li 

which upon optimizing in R yields 

llo^sTsILj < cilo^arsIlf/llsOoaTsIli? 

y y 

Upon taking L? norm in x, this proves (3.6). When combined with (3.2), we obtain from the above that 

\Kd^u\\L^ < C{t)^/^\\e^d^g\\]/^\\zeo.d^g\\^l 

y ^y ^y 

In order to prove (3.7), we use (3.3) with p = 1 and the ID Agmon inequality in the y variable to obtain 

< C{tfl^\\e^dTg\]^i^ {\\KdTdyg\\g, + ||8AST9IIl;)‘'" 

J^y ±Jy ±Jy Uy 

Taking the L? norm in x of the above yields (3.7), while an application of the ID Agmon inequality in x 
gives (3.8). For the v bounds, we use (3.3) with p = 2 and obtain 

l|8Tt.||i- < l|8T+‘«||i; < ||9„8r‘«lli~IIC'llL. 

<c{t)^/^\\e^8;'+'uU^ 

<C{tf>‘\\0^8T*'ghi. 

Integrating in x the above implies (3.9). An extra use of the ID Agmon inequality yields (3.10). □ 

Remark 3.2. The first two estimates in Lemma 3.1 also hold in the case when we don’t use the weight 
(i.e., when 9a = 1). Indeed, we use the relation 


ry 

u{y) = / giy) 
Jo 


exp 


y - y 
4(f) 


dy, 


which implies that 


for 1 < p < oo. 


Hy)\ < 
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3.2. Weighted Sobolev energy estimates for the good unknown. Let m > 0. We apply 5^ to (2.20), 
multiply the resulting equation with and integrate over H to obtain 

= -E(“) /ar'«‘«+‘Aar9-E(”) [ BT-’oeadyBige^a^g 

j=o \3 J J 

+ (”) / Blve„d^-ius„a^g 

= Um + Vm + Tm- (3-11) 


Here we have used the boundary conditions (2.21)-(2.22) and the cancellation 

J ^dT^'^gOld'^gdx = 0 , 

which follows upon integration by parts and the fact that dxi^fO"^) = 0. Dividing (3.11) by \\0ad^g\\L2, 
multiplying by r^Mm, and using the notations (2.26)-(2.30) and (2.35)-(2.37), we arrive at 


dt 


3“ Dr. 


, a(l-2a)~ 


+ 


2 — a 

'W 


— 


T^Mm 

Pad^gWL^ 


(Um + Kn + Tm) 


(3.12) 


In the next subsection we obtain lower bounds for the dissipative and damping terms on the left side of 
(3.12), while in the following subsection we estimate the nonlinear terms on the right side of (3.12). 


3.3. Bounds for the dissipative and damping terms. 

Lemma 3.3 (Poincare inequality with gaussian weights). Let g be such that dyg\y=Q = 0 and 
g\y=ggo = 0. Fov « E [l/4,1/2], m > 0, and t > 0 it holds that 

^^Pad^gWl^ < Po^dyd^gWl. (3.13) 


where 9a{t,y) = exp 


Proof of Lemma 3.3. The above inequality is classical, and it is a special case of the Treves inequal¬ 
ity which can be found in [Hor83]. For simplicity, we give a short proof for the case m = 0. Note that 

Oadyg = dy{9ag) - ^^^g 

as can be checked directly. Using that (a — 6)^ = (a + 6)^ — 4ab it then follows that 


J {9adygf dy = J (^dy{9ag) - dy 

= J (^dy{9ag) + ^dag^ dy - J 2y{9ag)dy{9ag) dy 
= J [dy{9ag) + dy + -^ J {9a,gf dy 

jdy, 

upon integrating by parts with respect to y in the third equality. No boundary terms arise in this process. □ 
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Using Lemma 3.3 we may bound the dissipation term in (3.11) from below as 




> 


PPadyd^gWl^ aV2/5 


a{l- (3) 




(3.14) 


where /3 G (0,1/2) is to be chosen precise later. 

For the damping terms in (3.11) we have the lower bounds 


1 


\\0ad^g\\L2 

1 


a(l — 2 a) 

4(t) 


2 — a, 


\WTg\\h + ^Il«„9™9iri2 


“ "•'1 - 2a)z^ + i;)'i^s„aTg\\h + ^ “(1 ‘"'' \Ka'^g\\h 


> 


\\0ad^g\\L^ 

a{l-2a)\\z9^dTg\\l2 , a^^/\l-2aY/^ 


+ 


(t) 

\\z9^d:^gh2 + 


8 (f) \\9^d^g\\L^ 4(f) " - """ (f) 

In the last inequality above we used that 

((1 - 20)^2 + 47 )^/^ > 2^/2 

((1 - 2 a ) z ‘^ + 47)^/^ > (1 _ 2a)^/22 

which holds for all z > 0, when a G [1/4,1/2] and 7 G [0,1/2]. In summary, in this subsection we have 
proven the following bounds. 

Lemma 3.4 (Lower bounds for the damping and dissipative terms). Fix a G [1/4,1/2], and let 
/3 ,7 G [0,1/2] be arbitrary. Then we have 


m>0 


(Dm + 


a(l - 2 a) 7 , 2 - a 

“r ^ 




4(f) ^ 2(f) 

a(l — 2 a),, ,, 


- ' 8 (f) 


+ 


a^/^/3 II II a7^/^(l — 2a)^/^ II „ 

H- —YJ^WdWDr.c H- —\ -Ilfl'IIZr.c 


2{ty 


4(f) 


l + a(l/2-7-/3)|| II 

+-(jj-Nix, 


(3.16) 


independently ofr > 0 . 

Proof. The lemma follows upon recasting (3.14) and (3.15) as 

b >^b 

2 2 (f)^/^ 


{t) 


2 — a a(l — 2 a) 

^ ™ + 4(f) 

and summing over m > 0 . 




a(l- 2 a)~ ^ ay^D(i- 2 ayD ^ ^ 

' I , i,\ “T 


8 (f) 


4(f) 


1 — a /2 — a 7 




□ 


3.4. Bounds for the nonlinear terms. In this subsection we bound the nonlinear terms on the right 
side of (3.12) for every m > 0, cf. estimates (3.20), (3.21), (3.22) below. When summed over m > 0 we 
obtain the following tangentially analytic estimates for the nonlinear terms. 
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Lemma 3.5 (Estimates for the nonlinearity). 

bounds 


There exits a universal constant C > 1 such that the 


E 

m>0 


\UjT^Mm 


c{tf^ . 

\\9\\X.r,A9\\Yr,c, 


E 

m>0 




^ Cjtf^ 

- r(i)V2 


\Y^. 


^ ^ Co{t)^^\ 1/2 

Pad^gh^ - r(t)V2 



lY-r.cc + 


£^!)!^!||„||V2 II ||1/2 


lY-r. 


(3.17) 

(3.18) 

(3.19) 


hold for every r > 0 and a € [1/4,1/2]. 

Proof. First, using (3.4)-(3.5), and the ID Agmon inequality in the x variable we obtain 



where C > 0 is independent of a E [1/4,1/2]. Upon multiplying by and using the definitions 

(2.26)-(2.30), the above bound implies 


C(t)V^ ('ft' yl/2yl/2 

Il«»8r9llz,» - (t(*))'/2 y >+■ J+2 

Similarly, by appealing to (3.9) and (3.10) we have 


+ E I ■ 

j = [m/2]+l 


(3.20) 




Padyg\\L2 


< 


[m/2] 

E 

j=o 


m 


jm-j 


v\\LlL<^Padydig\\L^i^2+ ^vh^ypadyd^xgh^. 


/=[m/2]+l 


[m/2] 




3=0 

m 2 ^ 

i=[m/2]+l ^ 


where C > 0 is independent of a E [1/4,1/2]. Upon multiplying by and using the definitions 

(2.26)-(2.30), the above bound implies 


\Vm\T^Mm ^ C{tf* 
Pad(fgh2 - (r(f))V2 





+ 


E 




1/2 


,y. 


1/2 


m—j+1^ m—j+2^3 


Di 


/=[m/2]+l 


(3.21) 
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For the last term on the right of (3.11) we appeal to (3.7), (3.8), (3.9), and (3.10) to obtain 


IKd^gh^ 


[m/2] , . 


+ 


2{t) 


E 


j=[m/2] + l 
mf2] 


m 


%v\\L2L^\\ead^ 


j^Q \J / 
[m/2] 




-i. 111/2 
L 2 

\ / / .y .y 

j =0 

m ^ V 

i=[m/2] + l ^ 

x||e„9™-^a,9ll[,/l|e„9™-"+‘8„9lll/ 

m /V 

+ Y. r)\\(^c.di-^"9\\L49ad:^-^9&^^^ 

i=[m/2] + l ^ 

where C > 0 is independent of a E [1/4,1/2]. Upon multiplying by r'^Mm and using the definitions 
(2.26)-(2.30), the above bound implies 




< 


\\9ad^9\\L^ - (r(f))i /2 


\ ^ y^i/2-^1/2 \x1/2 ^1/2 . \ ' -y vi/^ tri/^ ni/^ r)i/"i 

/ X j+1 ■*^j+2^m—j'^m—j / X ■*^i+l^m—j^m—j+l'^m—j'^m—j+1 

i=0 j=[m/2]+l 


+ 


C{t) 


1/4 /[™./2] 


(t(()) 1/2 


\ ^ '\^l/2-yxl/2 -^1/2 ^1/2 \ ^ wl/^ 1^1/^ yl/'l 7!/^ 

/ y ■*^/+2^m—j'^m—/ ”*” / X J+l^m—/^m—jr+l'^m—j'^m—/+! 

/=0 j=[m/2]+l 


(3.22) 


The proof of the lemma is completed upon summing (3.20)-(3.22) over m > 0 and using the bound 

m 

y~! Yh ^j^rn-j < Y1 “/ ^ 

m>0/=0 />0 A:>0 


for positive sequences {a/}j>o and {bj}j>o. □ 

Remark 3.6 (Analytic product estimates). We note that the proof of Lemma 3.5 directly implies that 
the following bounds hold 

I|u(9<‘')A9<"'IIx.,. < :;7|i7jll9''>llB„Jl9®lbv.. 

I|V(9'‘'H9®IU„. < :;^ll9®llfl„.ll9“>lliv,. 

4^I|V(9'")U(9'"')ILv„. < :;//72ll9'">llB,.Jl9'‘>ll59.. 
for some universal constant C > 0, independent of r > 0 and a E [1/4,1/2]. 
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3.5. Conclusion of the a priori estimates. At this stage we make a choice for the free parameters a, jS, 
and 7. First, we introduce 

6 = S{e) G (e, 1/10) 

which is to be chosen at the end of the proof, where without loss of generality e < 1/200. We set 


1-5 


a = 


13 = j = -. 
^ ' 2 


With this choice of a, /3, 7 , we sum estimate (3.12) for m > 0, appeal to Lemmas 3.4 and 3.5, and arrive at 

d „ „ 5/4 — 5 


llsll 




+ 


< + 






'Dr 


+ W^'IIsIId.,.) j llsllvv,. 


for some sufficiently large universal constants Cq, Ci > 1 which are independent of a and 5. Upon recalling 
the notations (2.34) and (2.37), we can rewrite the above in a more compact form as 


d II II 5/4 - 5 5 

■77 5 Wc “I-777— \\ 9 \\ Xr,a + 7777 

dt {t) Ci{t) 


5/4 


Br.a. + 


j ^ + 


^^0 


-(f)V 2 


\9\\Br 


\\Yr,c 

(3.23) 


with Co, Cl > 1 are universal constants, that are in particular independent of the choice of 5 G (e, 1/10). 
We next choose the function T{t) such that 


d 


:,T(r(t))=‘''"+3C„||s(t)l| 


dt 




T(t),Q: 


= 0 . 


(3.24) 


The above ODE is meant to hold a.e. in time, since the time derivative of the monotone decreasing absolutely 
continuous function (in fact Holder 1/2 continuous) is only guaranteed to exist almost everywhere. With 
this choice of r in (3.24), we infer from the a priori estimate (3.23) that 


A 

dt 


5 / 4 -( 5 | 


Ur.c + 


Clity 


{\\9\\Br,c, + \\9\\Br,c 


+ 


-(f)V 2 


Br,cl|y|IU,c ^ 


< 0 


which integrated on [ 0 , t] yields 


+ M-v 


Is, 


t(s),q: 


ds 


+ Co 


I 




< ll5o||x,g,„ < \\9o\\Xr^^l/2 < 

From (3.25) it immediately follows that 


(3.25) 


ft 


'0 




r(s),Q; “ — J 

which combined with (3.24) shows that we have the lower bound 


eCo 


r 


(f)3/2 > ro^/2 _ 


25 


(3.26) 
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for all t > 0, where C 2 = GCqCi is a universal constant that is independent of 6. From estimate (3.26) we 
see that the radius of tangential analyticity obeys 


r(t) > I 


on the time interval [0, Tg], where 


{Te)^ = 


3/2 


eC2 


(3.27) 


and we recall that 6 = 6{e) G (e, 1/10) is yet to be chosen. 

In order to see that the monotone decreasing analyticity radius is a Holder 1/2 continuous function of 
time, we may use the bound (2.38), integrate (3.24) from ti to ^2^ where 0 < < ^2 < 2/ are arbitrary, and 

use the estimate (3.25), to obtain 

Jti 

1/2 


, 1/2 

\Xr.o 


<6Co sup 

te[o,Te] ^ 

< 




rt2 

/ (i) 

J u 


-5\ 


3 

■^T,a j 


1/2 


J u 


dt\ 


1/2 


by using that 2(i — 1 < 0. 

To conclude the proof, we let 


6 = e log - 
e 

which is a permissible choice if s is sufficiently small. In that case, from (3.27) we obtain 

1 

1\ Elogi / 


% = 


log 

C 2 


— 1 = exp 


elog 


log 


^0 log 7 


^2 


- 1 . 


(3.28) 


(3.29) 


It is clear from (3.29) that as long as Tq'^^ log 2 > C 2 e^, which is ensured by (2.39), then we have that 

Te > exp 


elog 7 


for all 0 < e < 1/200, which concludes the proof of the a priori estimates. 


4. Uniqueness 

Assume G A/2 to,o with ||fl'o||x2To,o ^ e. Let^l^l and be two solutions to the system (2.20 )-(2.22) 
evolving from po^ with tangential radii of analyticity and respectively, which obey the bounds in 
Theorem 2.2. We fix 6 as given by (3.28). 

Also, define r(f) by 




r(0) = 


In view of the estimate (2.43) for ^1^1 and the lower bounds (2.41) for and r^^l, we have that 


22 < r(t) < 22 < 

8 - V 2 - 4 - 2 


(4.1) 


(4.2) 


for all t G [0, T^]. 


2 
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We consider the difference of solutions g = which obeys 


1 




- - {vdyg^^^ + v^^^dyg) + + v^'^-^u) 


(4.3) 


and has initial datum gQ = 0. Here we also denote u = = U(5') and v = =\{g). 

Using estimates for the nonlinear terms as in Remark 3.6, similarly to (3.23) we arrive at 


d 


5/4-'5„_ 


< 'r(i) + 


(t) 

:|l9<‘'(t)ll 


-(f)V2 ' 


X^(t) + 


r(t) 


Ci{t) 


2Co 




with Co,Ci > 1 being universal constants. Since r(t) < and the X^- q, norm is increasing in r, we 

obtain from (4.1) that 

+ < 0 . 

On the other hand, using (2.33), (4.2), and (2.42) we may bound 

4lls®(i)ll ■ ' 






Il9<"'(t)ll 


X 


< 


Combining the above two estimates with (4.4) we arrive at 

4ii9(*)iu„„+^44^ii9(*)iu.„,+ 


To conclude we note that by the definition of in (3.27) we have that 

5 _ elogi ^ 32£Co(t)^ ^ 2eCo{t)^ 


2eCo(f)^ 


Br{t)- 


Cl Cl 

holds. From (4.5) and (4.6) we obtain 


A 

dt 


W(t) + 


3/2 


5/4 — (5 


(i) 

which concludes the proof of uniqueness since gn = 0. 


-(/)3/2 


W(t) < 0 


(4.5) 


(4.6) 


5. Existence 


Throughout this section we fix a = 1/2 — <5, where <5 = elog 4. We assume the initial datum go obeys 
Ibo 11x2,0,1/2 < where the pair (tq, e) obeys (2.39). 

We first prove the existence of solutions gC) to a parabolic approximation of the Prandtl equations, with 
the term —vd^g present on the left side of (2.20). These solutions are shown to obey uniform in v bounds in 
^ n ^ for a sequence of tangential analyticity radii tC). These radii obey tC) > ro/2 

for all t G [0, T^] and are moreover uniformly equicontinuous on this time interval, where Tg is given by 
(3.27), i.e. 


{Tef 


3/2 


log 7 


iT* 


(5.1) 


for a sufficiently large universal constant iT*. Moreover, gC) and tC'> are shown to obey (3.24). 

With these uniform in bounds we then show that the tC) converge along a subsequence to an analyt¬ 
icity radius r(f) > ro/2 on [0, Tg], and along this subsequence, the gC) are shown to be a Cauchy sequence 
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in the topology induced by H L\BrQ^a- By the completeness of L°°{L?‘{6a{t,y)dydx)dt) the 

existence of solutions to Prandtl in the sense of Definition 2.1 is then completed. 


5.1. A dissipative approximation. For u > 0 we consider the nonlinear parabolic equation 


+ (n(") + + ^5^"^ 

o 

II 

1 

(5.2) 

dygyy=o = 9yy=oo 

(5.3) 

= Dig^^y := 0-iiy) f g^''\y)0i{y)dy 

Jo 


(5.4) 

v^^\y)=\{gy-.= - r dy^\my- 

Jo 


(5.5) 


Our goal is to construct solutions with corresponding tangential analyticity radii so that uniformly 
in u > 0 we have the estimate 


sup 

te[o,re] 








ds 


+ K, 


f 


(s) 


5/4-5 


(^) 




(0 




where Ff* > 0 is a sufficiently large universal constant, and the radii (t) obey the ODE 


(5.6) 


(5.7) 


For u > 0, estimate (5.6) and the ODE (5.7), correspond to (3.25) respectively (3.24) for the limiting 
Prandtl system u = 0. Although the system (5.2)-(5.5) is parabolic, we detail the construction of and 

since the first order ODE (5.7) has a nonlinear term which convergences only once the radius has 

been constructed already to satisfy this equation. The method of constructing g^'^'> and draws from ideas 
employed [KTVZll, IKZ12] for the hydrostatic Euler equations. 

At this stage it is convenient to introduce some notation. Eet > 1. Similarly to (2.26)-(2.37), for 
/i: H —)• R and r > 0 define the weighted Sobolev norms 


N 

\\h\\xN = ^m{h,T), 

m=0 

N 

IMd? = Dm{h,T), 

m=0 

N 

\\h\\z/y = ^rn{h,T), 

m=0 

N 

W^Wyj^ = 

m=l 

N 

WHbN = Bm{h,T), 

771=0 

We will use frequently that the bound 






Dm {hi t) 

E Zm{h, T ) 
Xm{h,T 



N 

~ 'y ^ Djn[h, t) , 

771=0 

N 

771=0 


Ymjh, Ty 

^lXm-l{h,Ty 


N 

Byn{h,T). 

777 =0 




holds independently of > 1 and r > 0. 


(5.8) 


(5.9) 
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5.2. A two-step Picard iteration for the dissipative system. We define 

S^''\t)ho = 

to be the solution to the initial value problem to the linear part of (5.2)-(5.5), namely 

= o (5.lO) 

{t) 

= 0 = 1^=00 (5.11) 

h^%=o = ho. (5.12) 

Solving (5.10)-(5.12) on H with the Neumann boundary condition (5.11) at y = 0 may be done using an 
even extension across y = 0 and solving the problem (5.10) on with vanishing boundary conditions as 
|y| —)• oo. As such, an explicit solution formula for (t) may be obtained, though it will not be essentially 
used here. We note that if ho obeys the boundary condition (5.11), the solutions S^'^\t)ho automatically lie 
in 2^2,1,/? for any /3 < 1 (cf. Definition 2.1). 

Next, we set up a two-step Picard iteration scheme. For n = 0,1 we let 

5(o>-)(t) = 5(h-)(t) = SM(f)yo 


while forn > 2 we define to be the mild solution (obtained by the Duhamel formula for the semigroup 
of the linear initial value problem 


_ ^Q2jn,u) ^ ^^Q^gin,.) + ^g{n,u) 

V"’"^|,=o = 0 = y("’")|,=oo 


(5.13) 

(5.14) 

(5.15) 


The pairing of g^'^ and g^^ 


in (5.13) is motivated by the bounds guaranteed by Remark 3.6. 


5.3. Sobolev bounds and convergence of the Picard iteration. Let N be an integer such that N > ^. 
For the remainder of this subsection we fix this value of N and we shall ignore the u and N indices for g 
and r. We claim that there exists jv > 0, to be chosen later, and a sequence of absolutely continuous 
monotone decreasing functions 

r(-) =4"); [0,r,,Ar] ^ 
with t(”)(0) = 7to/4 such that the bound 



< 2e (5.17) 


5ro Ttq 
4 ’ 4 


(5.16) 


holds for all n > 1, and some universal constant K > 1. 

We prove (5.17) inductively on n. For n = 1 this bound follows immediately from the assumption 
||fl'o||x 2 T^(j 1/2 — dissipativity of In order to prove the induction step we proceed as follows. 

Since Mm-i/2 < mM^ < 2Mm-i for all m > 1, and there are no boundary terms when integrating by 
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parts in x, for all m > 0 one may use Remark 3.6 to derive an estimate similar to (3.23) for the system 
(5.13)-(5.15) which is 

d y 5/4 5 

dr-' ""vw 




II I II II 

9 Wb^. . + IIS' IIb^ , 

t(") rir 


(*> ■ A'(t)V4 

+ ^ fils'”'" IIb», +lls‘"'"llg'' ,) Il9'”’llr» ,+jll9'”'llf-»« 

V tW / rr 4 \(n) 




tW 




+ 


K 


(n) 

K 


Il9'”’ll« 


(n) 




t(^) 


< f^^'> + 




1/2 

0 


T, 


lls' 


n—1) I 




r("-l) 




l9'”’lly» 






Y 


2^K 

N H- T-p^ 

(n-l) ,1/2 


(5.18) 


r("-l) 


'0 '0 
for all n > 2, where ft” is a sufficiently large universal constant (in particular 5, N,n, r-independent). In the 
second inequality in (5.18) we have used several times that cf. (5.16) we have 


-(n) 


max 


< max f < 2^/2. 


o<|j|<iv I j o<|i|<Af V5 

The main difficulty lies in obtaining a u—independent bound for the last term on the right side of (5.18). 
First we notice that since 

XmihiT') 


m 


and > 1 we may estimate 
2^+^K 


1/2 






B‘ 


(n-l) 




^ u (yjV+i(g("-i),r("-i))2 1!;!!/^||>-2)||2 


- 8"" 

< g lls*'”'~^^llv^+l + 


(n-l) uN " "^^(n-1) 

y:-» ' wv 


(n-2) I 


IS' 


(n-l) I 


(n-l) 


\rN 

T-(n-l) 


ly 


iV+1 

(n-l) 


^8lV+4^2^^)l/4||^(n-2)| 


\x^ IIS^”' ^^IIrJV 


(n-l) I 


ViV 


At this stage, forn > 1 we chose to solve the first order ODE 

12 ^ a: 


fin) + 


.1/2 


n(n-l) "^^(n-l)/ " ' 4 


7ro 


= 0, rH(0) = ^ 


(5.19) 


(5.20) 


The key point here is that by the induction step, the functions i) and i) are known, and due to the 
estimates (5.17) we have that 


/‘||9'"-"(»)IIbb „ +ll9'"-"(») 

Jo n("-l)(B) 


, , 2eK{ty 

Rjv ds < -r- < oo 

n("-l)(B) 0 


(5.21) 
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for all t G [OjTj^tv]- Thus the existence of an absolutely continuous solution to (5.20) is immediate. 
For n = 0 we may simply let (t) = 7ro/4. Moreover, from (5.21) we have that (5.16) holds at least on 
[0, T^^n], with Te^AT defined by 


3/2 


127V+i^2 

logi 


{Te,Nf- 


(5.22) 


In fact we will a-posteriori show that the time interval can be chosen independently of > l/v, as the 
factor 12^ is superfluous. For the moment however, the bound (5.22) is good enough since it is independent 
of n > 0 (recall that for now N is fixed). 

We now combine the bounds (5.18) and (5.19) with the choice for made in (5.20), integrate on 
[0, t], and use the induction assumption (via the bounds (5.17)) to obtain that 



< 2e (5.23) 


holds for all t G [0, T^^n]- In the second to last inequality above we have used that 

max{8^A:eTo/^4-^+2} < 1, 

1/2 

which holds if > 2 and Tq e is less than a small universal constant (which was assumed in (2.39)). This 
concludes the proof of the n-independent bounds (5.16)-(5.17) for the 

In order to show that the Picard approximation converges, we next show that the difference 

.(n) ^ (n) _ (n-1) 


contracts exponentially in a suitable weighted Sobolev space, of order N 
n > 1 define the decreasing function (t) by 






5tq/4 


5Tg/4 


= 0 , 


1 in X. For this purpose, for 


^(n)(o) = (5.24) 
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which by the uniform in n estimate (5.17) obeys 

5ro 


for 


We measure the difference by 

sup ] + 


[Oire^jv] 


(") (t) 






(")(s) 


4r" (7 


iN 


+ 


(s) 
K 


’t(")(s) 


ds 


.1/2 


(*)"■"-' (||9'”-‘*(s)IIb»-; + I|9‘"-"(s)IIb»-; ) I|9'">IL»-. ds. 

V "Sro/4 "5ro/4y r(»)(s) 


We claim that the sequence An contracts, and prove that 

^n-l + ^n-2 


An < 


(5.25) 


for all re > 2. In order to establish (5.25) we consider the equation obeyed by g^^^ 

+ .Lgin) 

= - U(p("-2))5,5(-2) - - Y{g(^-^'>)dyg^^-^'> 

+ ^V(5("-i))U(5("-^)) + ^V(5("-^))U(p("-2)). 


Using estimates that are similar to those in Remark 3.6, the bounds (5.18)-(5.23), and using the choice of 
in (5.24) it then follows that 

, j_:_ii;^v'oii , j_I ii;rvvii ,, , j_ ii , i j_ ^' 




fll9<”'llB»-. + ll9'”'llB»-.)+jll9‘">ll 

V •f(") / 4 


\vN 

AM 


+ 


8^K 


'B, 


i-i + 
'5Tg/4 


3JV-1 I ll^^’^^llyiv-i 
ST-o/ry fM 


< 


dt 


.fin) + 


iN 


K 


U/2 


Is 


Sto/4 


6to/4 , 


Il9<”'lly»-. 


(n) 


+ 


1/2 I 




:(n-2) 


(ti-2) 16 


lyre 


+ 


A^K 
1/2 I 


+8«/U{«)‘/*||yn-2)|| ||yn-2)|| 


5to/4 


5to/4 


5tq/4 


< 


/V2^+^7t:, 


r, 


3/2 

0 


^^llnre-i ^^llyre ^ + — 


p("-2) 


(n-2) 16 


yre 


1 / 8^K 




5to/4 


B. 




ly- 


(n-l) 


ly- 


(ra-l) 


(5.26) 


for t G [0, Tg jv]- The proof of (5.25) now follows from (5.26) upon integrating in time, recalling that 
6 = clog i, that (e, tq) obey (2.39), and that the bound 


N2^+^K 

3J2 


Il9'”-">ll 




AeN2^K 

(n-2) “ 3/2, A 5/4-5 


(5 16eiV2^iT2(T<,,jv)'^ 




5/4 


3/2 X 
-l-Q ^ 


< 


4K{t) 


5/4 
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holds in view of the bound (5.17) and the definition of jv in (5.22). Thus, we have proven (5.25), from 
which it follows that 


0 < < ao 


VTt- l' 


+ 0,\ 


'V^+i 


0 as n 


oo 


where ao,ai > 0 are determined from computing Ai and A 2 . This concludes the proof of convergence 
for the Picard iteration scheme (5.13)-(5.14) on [OjTg^Tv]. The convergence holds in the norm defined by 
An- Moreover, fhe available bounds are sufficienf in order fo show fhaf fhe limiting funcfion obeys 
(5.2)-(5.5) poinfwise in x when infegrafed againsf H^{6ady) functions of y. 


5.4. A posteriori estimates for the dissipative approximation. Having constructed solutions of 
(5.2)-(5.5) with finite Sobolev regularity in x (of order N > 1 /v), we a posteriori show that these solutions 
obey better bounds, and in particular, are real-analytic with respect to x. 

For this purpose, we would like to perform estimates similar to those in the previous subsection, and 
pass N ^ 00 . The main obstruction to directly using the bound (5.17) and passing —)■ 00 is that the 
time of existence we have so far guaranteed for g^'^^ is defined by (5.22), and thus depends on N itself. 
Thus, the first step is to show that obeys i/-independent Sobolev bounds on a time interval Tg that is 
independent of N (and v). 

As before, let N be such that > 1 with the caveat that we will in this subsection look for bounds 
independent of N. Let K be the constant from (5.18), and define T^\t) by 



+ 119 M 

\^) 

N 


\b^, . 

t(^) 


= 0 , 


(5.27) 


with initial value t^^(O) = 7ro/4. For each N, this is a first order ODE, with a degree N polynomial 
nonlinearity in Due to the a priori bounds (5.17) inherited by g^'^\ at least on [0, T^^n] the ODE (5.27) 
has an absolutely continuous solution. We let be the maximal time for which stays above 5ro/4. 
On [0, T*j^] all the estimates in the previous section are justified. We already have shown that T* > T^^n, 
and we now claim that for the > 0 defined in (5.1), which is independent of N and u. 

With as defined above, and N > arbitrary, we perform an estimate in the spirit of (5.18)-(5.19), 
use the definition of in (5.27), and arrive at 
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where > 1 is a universal constant. Integrating the above on [0, T] and using that ||5 o||x 2 to — 
N > 1/u we obtain 


sup ( ^\\9^''Ht)\\xN 

te[o,T] \ 




r- 

lo (s)^ 


+ 


K r 

Jo 




-s 


(^) 


is 


, +11/ 

‘'l(s)|lBiv ) ds 

7(s) 


0 B'v 

+ II9("H^)IIb'v 

Phs) 



Ml 




ds 


< e. 


(5.29) 


Estimate (5.29) above implies that 
AK P 


1/2 


M 


rP(s) 


ds < 




6ti 


1/2 


^0^^ log I 


16K^ / (s)^/^||p('^l(s) 


/o 


lx' 




Ils'^’Wlh 


, 16it'3e2 iQK^e AK’^Uy 

ds < ---=-^ < 


(s) 




(s) 


log J log i 


upon appealing to (2.39). Inserted in (5.27), the above bounds a posteriori show that 


_(.) 7to 

Tx {t) ^ — -JJ^ ^ 

^ ^0 log 7 


> 


5ro 


for all f < Tg, as long as obeys 


(5.30) 


(5.31) 


3/2 


{Tef < 


log 


16X2 


It is clear that the Tg defined earlier in (5.1) obeys the above estimate if X* is taken sufficiently large. This 
shows that T* ^ for each N > Moreover, the bound (5.31) which combined with (5.29) yields 



< e, 


(5.32) 


for any X > 1, where X > 1 is a fixed universal constant. Note that upon passing X —oo in (5.32), and 
using the Monotone Convergence Theorem, we also obtain the bound 



< 4c, (5.33) 


for the real-analytic norms of <7(^1. Due to the monotonicity of the norms with respect to r, this proves (5.6). 

In order to obtain a limiting analyticity radius in the limit as X —)• 00, which obeys the nonlinear 
ODE (5.7), we may first try to show that the sequence of absolutely continuous functions 1^ 

in fact equicontinuous on the time interval [0, T^]. This seems however not possible due to the third term 
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on the left side of (5.27). We instead define a new sequence of radii for which the trick used to prove 
uniqueness in Section 4 applies, and we are able to prove that uniformly equicontinuous (in 

fact uniformly Holder-1/2 in time). Let 




-f 


2K 










= 0, 05^^(O)=to 


(5.34) 


The existence of solutions to (5.34) is immediate since the nonlinearity is a polynomial of finite degree, with 
coefficients that are integrable in time by (5.33). We next observe that in view of (5.30), by using a version 
of (5.31), we arrive at 


for all tG[0,T,]. 


Now, similarly to (5.28) we have that 
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By 1 


7^) 


(^)l 






By 1 


(5.35) 


In the last inequality of (5.35) we have used the same trick as in Section 4: that by (5.32) we have 


sup 

te[o,rd V Bo 


5to/4 


2>2K^e^ , 1 , 

<-< e log -= 8 

To e 


upon appealing to assumption (2.39). 

Using the bound (5.35), we show that the sequence of absolutely continuous functions {0^^}Ar>j^-i, 
is in fact uniformly bounded in C'^/^([0,Tg]), and thus uniformly equicontinuous. For this purpose, let 
ti,t 2 G [0, Tg] be such that |fi — 12 \ < C- Using the mean value theorem, the definition of in (5.34), and 
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the definitions (5.8)-(5.9) we arrive at 


Tq Jtl 
4:K /■*= 




(s)|| 5 iv ds 
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.1/2 


pt2 

/ (^) 
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l/8|lnM 
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(5.36) 


Since C £ (0,1) was arbitrary, it follows from (5.36) that the are uniformly equicontinous. The Arzela- 
Ascoli theorem guarantees the existence of a subsequence with ^ oo as k ^ oo, and of a function 
such that 

^ uniformly on [OjTg] as k^oo. 

Moreover, we have that > ro/2 on [0, Tg]. By passing N = Nk ^ oo in (5.34) we obtain 


dt 


+ 


II (y)|| 
1/2 11^ 


= 0 , 


r('")(0) = tq. 


(5.37) 


In order to justify (5.37) we use that by (5.33) we have that IIb 5 to /4 ^ ^e])> that the con¬ 
vergence of is uniform. Moreover, using (5.33) and a bound similar to (5.36), it follows from 

(5.37) that 





16Ke^ \ 
) 




uniformly for ti,t 2 G [0, Tg]. That is, the radii are uniformly (with respect to u) Holder 1/2 continuous. 
This concludes the proof of (5.6)-(5.7). 


5.5. Existence of solutions to the Prandtl system. It remains to pass u —)• 0 and obtain a limiting 
solution g of the Prandtl equations (2.20)-(2.22), in the sense of Definition 2.1, of a tangential analyticity 
radius r which solves (3.24), such that the pair (p, r) obeys the bounds (2.41)-(2.44). 

We have shown in the previous subsection that the sequence is uniformly equicontinuous, and thus 
by the Arzela-Ascoli theorem we know that along a subsequence —)• 0, we have that —)• r uniformly 
on [0, Tg], with r(0) = tq , and T{t) > tq/2 on this interval. Note that from the bound (5.33) it follows that 
(5.6) holds with replaced by the limiting function r. 

Without loss of generality, the above subsequence {i'k}k>i obeys 


Ufc+i i/fc A;2 ■ 

We next show that the subsequence is Cauchy in the norm induced by the left side of (5.6). For this 
purpose, let 

gk = - 5^'"'=+!^ 
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and define 


where 


Gk ■■= sup ^\\gk{t)\\x. 

lO,Te] ^ 


2K 


■||5/c(^)|| 


I / \ (5 

'o (s) 


ds 
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5/4-<5 


0 Tk 


(S)V2 
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Stq/4 


Y- , \ds 


d _ ‘2K 
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'5to/4 


= 0, rfc(O) = 


6ro 


(5.38) 


Note that from the bound (5.33), upon choosing sufficiently large, we obtain that 

fk{t) > To for all t G [0, T^]. 

We claim that 

1 


Gk < 




(5.39) 


To prove (5.39), we consider the equation obeyed by gk 

dm - dfiQk + Kipd^gk + -^gk - T^kdlSk 

= {^k - k'k+i)dlg^''^+d - \]{g^''>^"i)d^gk - ^{gk)dxg^''’^+^'’ 


- - Y{gk)dyg^^^'^ + {gk) + ^V{gk)Y{g^^^^d)_ 


Similarly to (4.4) we obtain that 
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(5.40) 


In the last inequality above we used that (5.38) implies that Tk < 6ro/5 < 5ro/4, and thus 


sup 

, m>0 
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V5ro/4 
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upon possibly increasing the value of K. Using the definition of we obtain from (5.40) that 
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At this stage, we use (5.33) and (5.1) which imply that 

/ 2K^ 


sup 

te[o,re; 


6ti 


3/2 




X, 


< = 




5^0/4 1-3/2 

/ <5^0 


.3/2 


log 


-{Ter < 


1 


'0 / "'0 '0 

if is sufficiently large. Then, again appealing to (5.33), upon integrating (5.41) in time we obtain 

Gk<^< ^ 


fc2rn k?Tn 


From here it follows that is a Cauchy sequence in the topology induced by G^- 

Thus, there exists a limiting function g € L°°([0, T^]-,Xt^) n T^([0, Tg]; Br^ n Bt-q) such that g 

in this norm. In particular, it immediately follows that g G L°°([0, T'e]; ^2,1,Q:/{t>) thus 5 is a solution of 
the Prandtl equations (2.20)-(2.22) in the sense of Definition 2.1. 

Finally, using these bounds one may show that upon passing u = —>• 0 in (5.37), the limiting 

analyticity radius r(f) and the solution g{t) obey the ODE (3.24). This concludes the proof of the existence 
of solutions to (2.20)-(2.22). 


Acknowledgments. The authors are thankful to Igor Kukavica, Toan Nguyen, and Marius Paicu for 
stimulating discussions, and to the anonymous referees for helpful suggestions. The work of VV was sup¬ 
ported in part by the NSF grants DMS-1348193, DMS-1514771, and an Alfred P. Sloan Fellowship. 

Conflict of Interest. The authors declare that they have no conflict of interest. 


[AWXY14] 

[CGPll] 

[CKV14] 

[CLSOl] 

[CSOO] 

[DR04] 

[EE97] 

[GGN14a] 

[GGN14b] 

[GGN14c] 

[GNU] 

[GN14] 

[GreOOa] 

[GreOOb] 

[GSS09] 

[GVDIO] 

[GVM13] 

[GVN12] 


References 

R. Alexandre, Y.-G. Wang, C.-J. Xu, and T. Yang. Well-posedness of the Prandtl equation in Sobolev spaces. J. Amer. 
Math. Soc., 2014. 

J.-Y. Chemin, I. Gallagher, and M. Paicu. Global regularity for some classes of large solutions to the Navier-Stokes 
equations. Ann. of Math. (2), 173(2):983-1012, 2011. 

P. Constantin, I. Kukavica, and V. Vicol. On the inviscid limit of the Navier-Stokes equations. arXiv:1403.5748, Proc. 
Amer. Math. Soc., 2014. 

M. Cannone, M.C. Lombardo, and M. Sammartino. Existence and uniqueness for the Prandtl equations. C. R. Acad. 
Sci. Paris Ser I Math., 332(3):277-282, 2001. 

R.E. Caflisch and M. Sammartino. Existence and singularities for the Prandtl boundary layer equations. ZAMM Z. 
Angew. Math. Mech., 80(ll-12):733-744, 2000. 

P.G. Drazin and W.H. Reid. Hydrodynamic stability. Cambridge University Press, 2004. 

W. E and B. Engquist. Blowup of solutions of the unsteady Prandtl’s equation. Comm. PureAppl. Math., 50(12): 1287- 
1293, 1997. 

E. Grenier, Y. Guo, and T. Nguyen. Spectral instability of characteristic boundary layer flows. arXiv:1406.3862, 2014. 
E. Grenier, Y. Guo, and T. Nguyen. Spectral instability of symmetric shear flows in a two-dimensional channel. 
arXiv:1402.1395, 2014. 

E. Grenier, Y. Guo, and T. Nguyen. Spectral stability of Prandtl boundary layers: an overview. arXiv:1406.4452, 2014. 
Y. Guo and T. Nguyen. A note on Prandtl boundary layers. Comm. Pure Appl. Math., 64(10): 1416-1438, 2011. 

Y. Guo and T. Nguyen. Prandtl boundary layer expansions of steady Navier-Stokes flows over a moving plate. 
arXiv:1411.6984, 2014. 

E. Grenier. On the nonlinear instability of Euler and Prandtl equations. Comm. Pure Appl. Math., 53(9): 1067-1091, 

2000 . 

E. Grenier. On the stability of boundary layers of incompressible Euler equations. J. Differential Equations, 
164(1): 180-222, 2000. 

E. Gargano, M. Sammartino, and V. Sciacca. Singularity formation for Prandtl’s equations. Phys. D, 238(19): 1975- 
1991, 2009. 

D. Gerard-Varet and E. Dormy. On the ill-posedness of the Prandtl equation. J. Amer. Math. Soc., 23(2):591-609, 

2010 . 

D. Gerard-Varet and N. Masmoudi. Well-posedness for the Prandtl system without analyticity or monotonicity. 
arXiv:1305.0221, 2013. 

D. Gerard-Varet and T. Nguyen. Remarks on the ill-posedness of the Prandtl equation. Asymptotic Analysis, 77:71-88, 

2012 . 







30 M. IGNATOVA AND V. VICOL 

[HH03] L. Hong and J.K. Hunter. Singularity formation and instability in the unsteady inviscid and viscous Prandtl equations. 
Commun. Math. Sci., 1(2):293-316, 2003. 

[Hor83] L. Hormander. The analysis of linear partial differential operators III, volume 257. Springer, 1983. 

[IKZ12] M. Ignatova, I. Kukavica, and M. Ziane. Local existence of solutions to the free boundary value problem for the 
primitive equations of the ocean. Journal of Mathematical Physics, 53:103101, 2012. 

[Kla83] S. Klainerman. On “almost global” solutions to quasilinear wave equations in three space dimensions. Comm. Pure 
Appl. Math., 36(3):325-344, 1983. 

[KMVW14] I. Kukavica, N. Masmoudi, V. Vicol, and T.K. Wong. On the local well-posedness of the Prandtl and the hydrostatic 
Euler equations with multiple monotonicity regions. SIAM J. Math. Anal, 46(6):3865-3890, 2014. 

[KTVZll] I. Kukavica, R. Temam, V. Vicol, and M. Ziane. Local existence and uniqueness for the hydrostatic Euler equations 
on a bounded domain. J. Differential Equations, 250(3): 1719-1746, 2011. 

[KVll] 1. Kukavica and V.C. Vicol. The domain of analyticity of solutions to the three-dimensional Euler equations in a half 
space. Discrete Contin. Dyn. Syst., 29(l):285-303, 2011. 

[KV13] 1. Kukavica and V. Vicol. On the local existence of analytic solutions to the Prandtl boundary layer equations. Com¬ 

mun. Math. Sci., ll(l):269-292, 2013. 

[LCS03] M.C. Lombardo, M. Cannone, and M. Sammartino. Well-posedness of the boundary layer equations. SIAM J. Math. 
Anal, 35(4):987-1004 (electronic), 2003. 

[LWX15] W. Li, D. Wu, and C.-J. Xu. Gevrey class smoothing effect for the Prandtl equation. arXiv:1502.03569, 02 2015. 

[LWY14] C.-J. Liu, Y.-G. Wang, and T. Yang. A well-posedness theory for the Prandtl equations in three space variables. 

arXiv:1405.5308, 05 2014. 

[Mae 14] Y. Maekawa. On the inviscid limit problem of the vorticity equations for viscous incompressible flows in the half¬ 
plane. Comm. Pure Appl. Math., 67(7): 1045-1128, 2014. 

[MW12] N. Masmoudi and T.K. Wong. On the theory of hydrostatic Euler equations. Arch. Ration. Mech. Anal, 
204(1):231-271, 2012. 

[MW14] N. Masmoudi and T.K. Wong. Local-in-time existence and uniqueness of solutions to the Prandtl equations by energy 
methods. arXiv:1206.3629, Comm. Pure Appl. Math.,, 2014. 

[01e66] O.A. Oleinik. On the mathematical theory of boundary layer for an unsteady flow of incompressible fluid. J. Appl. 
Math. Mech., 30:951-974 (1967), 1966. 

[OTOl] M. Oliver and E.S. Titi. On the domain of analyticity of solutions of second order analytic nonlinear differential 
equations. J. Differential Equations, 174(l):55-74, 2001. 

[Pra04] L. Prandtl. Uber Elussigkeitsbewegung bei sehr kleiner Reibung. Verh. Ill Intern. Math. Kongr. Heidelberg, Teuber, 
Leipzig, pages 485^91, 1904. 

[PVll] M. Paicu and V. Vicol. Analyticity and gevrey-class regularity for the second-grade fluid equations. J. Math. Eluid 
Mech., 13(4):533-555, 2011. 

[PZ14] M. Paicu and Z. Zhang. Global well-posedness for 3D Navier-Stokes equations with ill-prepared initial data. J. Inst. 
Math. Jussieu, 13(2):395^1L 2014. 

[SC98] M. Sammartino and R.E. Caflisch. Zero viscosity limit for analytic solutions, of the Navier-Stokes equation on a 
half-space. 1. Existence for Euler and Prandtl equations. Comm. Math. Phys., 192(2):433^61, 1998. 

[WXY14] Y.-G. Wang, E. Xie, and T. Yang. Local well-posedness of Prandtl equations for compressible flow in two space 

variables. arXiv:I407.3637, 07 2014. 

[XZ04] Z. Xin and L. Zhang. On the global existence of solutions to the PrandtPs system. Adv. Math., 181(1):88-133, 2004. 

[ZZ14] P. Zhang and Z. Zhang. Long time well-posdness of Prandtl system with small and analytic initial data. 

arXiv:1409.1648, 2014. 

Department of Mathematics, Princeton University, Princeton, NJ 08544 
E-mail address: ignatovaOmath. pr inceton . edu 

Department of Mathematics, Princeton University, Princeton, NJ 08544 
E-mail address: vvicol@math. princeton . edu 



